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ABSTRACT 

The formation process of circumstellar disks is still controversial because of the interplay of complex 
physical processes that occurs during the gravitational collapse of prestellar cores. In this study, we 
investigate the effect of the Hall current term on the formation of the circumstellar disk using three- 
dimensional simulations. In our simulations, all non-ideal effects as well as the radiation transfer are 
considered. The size of the disk is significantly affected by a simple difference in the inherent properties 
of the prestellar core, namely whether the rotation vector and the magnetic field are parallel or anti¬ 
parallel. In the former case, only a very small disk (< I AU) is formed. On the other hand, in the 
latter case, a massive and large (> 20 AU) disk is formed in the early phase of protostar formation. 
Since the parallel and anti-parallel properties do not readily change, we expect that the parallel and 
anti-parallel properties are also important in the subsequent disk evolution and the difference between 
the two cases is maintained or enhanced. This result suggests that the disk size distribution of the 
Class 0 young stellar objects is bimodal. Thus, the disk evolution can be categorized into two cases 
and we may call the parallel and anti-parallel systems as Ortho-disk and Para-disk, respectively. We 
also show that the anti-rotating envelopes against the disk-rotation appear with a size of > 200 AU. 
We predict that the anti-rotating envelope will be found in the future observations. 


1. INTRODUCTION 

Circumstellar disks are born around protostars in the 
course of the self-gravitational collapse of the molecular 
cloud cores. The angular momentum evolution during 
the collapse is critical for the disk formation because the 
disks are supported by centrifugal force. 

The magnetic field plays a central role for the angular 
momentum evolution. During the gravitational collapse, 
a toroidal magnetic field is created by the rotation, and 
the magnetic tension decelerates the gas rotation, remov¬ 
ing the angula r momentum. This effect is known a s mag¬ 
netic braking ([Mouschovias and PaleologoullI979| l. Pre¬ 
vious studies using ideal magnetohydrodynamics (MHD) 
simulations have shown that, for a typical magnetic 
field strength, the disk formation is completely sup¬ 
pressed in the early Class 0 phase (iMellon and Lil 120081 
iHennebelle and Fromandr2008l iBate et a/.! 120141 1. 

Meanwhile, the observations of Class 0 young stellar 
objects (YSOs) showed that a relatively large circum- 
stellar disk (r ^ 50 AU) exists around young protostar s 

(|Ohashi et aLll20ll : iSakai et '^12014 iTobin et aLll2015ll . 

These observations indicate that some very young proto¬ 
stars have relatively large circumstellar disks with a size 
of r > 10 AU, suggesting a disagreement between the 
previous theoretical works and the observations. 

Possible physical mechanisms for resolving the discrep¬ 
ancy between theory and observation are the Ohmic and 
ambipolar diffusion (IMel lon and Lill2009t iMachid a et al\ 

I20l4l^mida et ^l2015HTsukamoto et aLll2015al l. Pre- 
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vious works following the formation of protostars showed 
that a small disk with size of r < 1 AU is formed around 
the protostar when the magnetic diffusion are considered 
(|Tomida et aLl[2M^ iTsukamoto et^l2015ah . However, 
the formation of a disk with a size of r > 10 AU at the 
early phase of protostar formation in a typically magne¬ 
tized cloud core is still highly difficult even with these 
effects. 

The effect of the Hall current term is the least stud¬ 
ied effect in the context of disk formation. The Hall 
current term generates a toroidal magnetic field from 
a poloidal magnetic field and directly affects the mag¬ 
netic tension that determines the magnetic braking effi¬ 
ciency. The Hall current term in MHD equations is not 
invari ant against the globa l inversion of the magnetic 
field (iWardle and Nd 119991: iKrasnopolsky et aLl 120111 : 
iLi et al\ 120111 : [Braiding and Wardlel I2012D and its ef¬ 
fect changes depending on whether rotation vector and 
magnetic field of host cloud core are parallel or anti¬ 
parallel. When bf the Hall diffusion coefficient is neg¬ 
ative (this is true for p < 10“^^ g cm“^ in our model, 
see figure |4|) and the rotation vector and the magnetic 
field are anti-parallel, the Hall current term weakens 
the magnetic braking. Meanwhile, the Hall current 
term strengthens the magnetic braking in the parallel 
case. Despite the possible importance of the Hall cur¬ 
rent term in the disk evolution, it is still unclear how 
the Hall current term affects the formation of the cir- 
cumstellar disk because the previous nu merical studies 
(|KrasnoDolskv et ^ 120111 : iLi et al\ 1201111 neglected the 
first core evolution phas e, which plays an important role 
in the disk formation (I Mac hida and MatsumotcJ 120111: 
Dapp^et ^l2012t iTomida ~t oi.! 120151 : iTsu^moto et oil 
2015ajl and simplified the radiative transfer. Three- 
dimensional simulations are also necessary for investi¬ 
gating non-axisymmetric effects. 

In this paper, we performed three-dimensional simula¬ 
tions starting from prestellar cloud cores. Our numerical 
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simulations include all non-ideal MHD effects as well as 
the radiative transfer. The simulations were conducted 
until the birth of the protostar. We did not use any sink 
technique for the center and hence our simulations do not 
suffer from numerical artifacts introduced by sink parti¬ 
cle or the inner boundary which may art ificially change 
the f ormation and evolution of the disk (|Machida et al\ 

[Ml . 

2. NUMERICAL METHOD AND INITIAL CONDITIONS 


In this study, we solved the non-ideal radiation 
magneto-hydrodynamics equations with self-gravity. 
The numerical method, except for the Hall cur¬ 
rent term, is the same as t h at use d in our pre¬ 
vious study ([Tsukamoto et al\ l2015a|l . The ideal 
MH D part was solved u si ng th e met hods proposed 
by llwasaki and Inutsukal (120111 I2013I1 . The ra- 
diative transfer was treated with the methods o f 
([Whitehouse and Batel 12004 iWhitehouse et al\ 1200511 . 
We treated the Ohmic and ambipolar diffus io n with 
the method described b y iTsukamoto et al\ (l2013a[l 
and IWurster et al\ (|20l4l . respectively. Both diffu- 
sion processes were ac celerated by super time stepping 
([Alexiades et a/.lll996ll . 

For this study, we n ewly implemen t ed th e Hall cur¬ 
rent term according to IWurster et al\ (120141 . The Hall 
current term couples with the ideal terms and changes 
the phase velocity. Thus, it is unclear that the sub- 
cycle method, which is often used for the diffusi on terms 
([Machida et aZ.ll^ll4 ITsukamoto et aLll2015all . is valid 
for the Hall current term. Therefore, in our simula¬ 
tions, the ideal MHD term and the Hall current term 
are updated with the same time-step. The smaller time- 
step of either the ideal te rm or the Hall curren t term, 
^%aii = Chaii^^/(47r|?7H|) (ISano and Stonell2002ll is used 
for the update. Here, Chaii = 0.4 is the Courant- 
Friedrichs-Levy number for the Hall term, and h is the 
smoothing length. 

We conducted the numerical tests for the Hall current 
term and confirmed that the scheme can correctly cal¬ 
culate the whistler mode in the linear wave propagation 
test. Furthermore, we conducted gravitational collapse 
tests of non-rotating cloud cores and confirmed that the 
rotation amplitude induced by the Hall current term does 
not depend on the direction of the magnetic field. We 
adopted the equati on of state (EOS), a n d dust and gas 
opacit y tab les fromiTomida et al. (1201311 . ISemenov et al\ 
([200311 . and iFerguson et aTT i 2005 1. respectively. We em- 
ployed the resis t ivity table used in our p revious works 
([Okuzuniil 120091 : ITsukamoto et al\ l2015al) in which the 
fixed dust grain size of a = 0.035yLtm and fixed cosmic- 
ray ionization rate of ^CR = are assumed. 

We modeled the initial cloud core with an isother¬ 
mal uniform gas sphere using about 3 x 10® particles. 
The mass and temperature of the initial core are 1 
Mq and 10 K, respectively. Initially, the core has a 
radius oi R = 3.0 x 10® AU and is rigidly rotating 
with an angular velocity of Dq = 2.2 x 10“^® s“^. 
The initial magnetic field is uniform and parallel to 
the rotation (z-) axis with a magnitude of Hq = 1.7 x 
10^/rG. The corresponding initial mass-to-flux ratio rel¬ 
ative to the critical value is ^ = (M/$)/(M/$)crit = 4 
where $ = and = (0.53/37r)(5/G)^/® 

([Mouschovias and Spitzenl 197611 . 


We conducted three simulations. Model Para, Ortho, 
and NoHall. The magnetic field and the rotation vector 
are initially perfectly parallel in Model Ortho, NoHall 
and anti-parallel in Model Para. The model NoHall does 
not include the Hall current term. Other parameters are 
the same in the models. The runtime of Models Para. 
Ortho, and NoHall was about 3.7x 10®, 1.4 x 10®, 1.6 x 10"^ 
CPU hours, respectively with XC30 in NAOJ. 

A boundary condition is imposed at i?out = 0.995i?, 
and the particles with r > i?out rotate with an initial an¬ 
gular velocity. Thus, the gas is confined in a rigidly ro- 
tating shell. This boundary is ver y similar to that used in 
iMatsumoto and Tomisakaf ([200411 : iMachida et al\ (1200711 
and a lso used in our previous work (ITsukamoto et al\ 
l2015all . In addition, a boundary condition for radiative 
transfer is introduced by fixing the gas temperature to 
be 10 K when p < 2.0 x 10”^^ g cm“®. 


3. RESULTS 

In figure (U we show the structure at the center of the 
cloud core. The left and right panels show the result of 
Models Para and Ortho, respectively. The central den¬ 
sities are 10“® g cm“® for the left panels, which cor¬ 
respond to slightly before the protostar formation, and 
10“^ g cm“® for the right panels, which is immediately 
after the protostar formation. 

When the rotation vector and magnetic field are in 
the anti-parallel configuration, a large disk with a size 
of r ~ 20 AU is formed (Model Para; top left panel). 
The disk is so massive that spiral arms are created by 
the gravitational instability. We confirmed that the 
Toomre’s Q value, Q = /«epCs/(7rGE) is Q ~ 1 in 
the entire disk region (5 < r < 20 AU). In the top 
panel of figure O we show the force balance between 
the pressure gradient force, the centrifugal force, and 
the radial gravitational force on the x-axis. The gas 
is mainly supported by the centrifugal force. There¬ 
fore, a rotationally-supported massive disk is formed in 
Model Para. On the other hand, when the rotation vec¬ 
tor and magnetic field are in the parallel configuration, 
no large disk (r > 10 AU) appears (Model Ortho; top 
right panel) because the magnetic braking is strength¬ 
ened by the Hall current term and the angular momen¬ 
tum is efficiently removed from the central region. The 
dense region {p > 10“^^ g cm“®) at the center, which 
has a radius of r ^ 5 AU, is the remnant of the first 
core and is not a rotationally-supported disk. Although 
a rotationally-supported disk with r < 0.6 AU is also 
formed in Model Ortho around the protostar, as shown 
in the bottom panel of the figure H] the difference of the 
disk size of Model Para and Model Ortho is remarkable. 

The plasma (3 in the disk regions (r < 20 AU) of Model 
Para is large {j3 > 100) because the magnetic flux is 
largely removed due to the magnetic diffusion. In this 
high /3 region, the magnetic field and the gas are almost 
decoupled and the magnetic braking is no longer impor¬ 
tant. The disk is sufficiently massive and develops grav¬ 
itational instability. Thus, gravitational instability may 
play an important role for angular momentum transfer 
in the subsequent evolutionary phase. Furthermore, in 
such a massive extended disk, disk fragmentation, which 
is a promising mech anism for the format ion of binaries or 
wide-orbit planets, ([Marois et possibly occurs 

in the subsequent evolution ( BossI 11994 llnutsuka et ai\ 
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Fig. 1.— Cross sections of the density (top panels) and the plasma /3 (bottom panels) in the x-y plane. The left panels show the results of 
Model Para, while the right panels show those of Model Ortho. The central densities at these snapshots are 10“® g cm~^ and 10“^ g cm“^ 
for the left and right panels, respectively. 


20K; iMachid a et al\ l2011bt iTsukam oto and Machidal 
20nl l201,lHTsukamoto erd\\2mWiJ2()lF>U\ . Thus, the 
parallel or anti-parallel property of the cloud core would 
play a crucial role for the formation of the binary or 
wide-orbit planets. 

To quantify the strength of the rotation at the center 
of the cloud core, we show the mean specific angular mo¬ 
mentum of regions with p > g cm“^ as a function 

of the central density in figure [H The figure shows that 
the specific angular momentum in Model Para is about 
an order of magnitude larger than that of Model Ortho. 
The specific angular momentum in Model Ortho (Para) 
is about three times smaller (larger) than that in Model 
NoHall. The combination of the spin-up effect (weak¬ 
ening of the magnetic braking) in the anti-parallel case 
and the spin-down effect (strengthening of the magnetic 
braking) in the parallel case causes the large difference. 

The mass and the magnetic flux $ = J BdS of re¬ 
gions with p > g cm“^ in model Para, Or¬ 

tho, and NoHall at the beginning of the second col¬ 
lapse {pc = 10“® g cm“^), was (M (Mq), $ (G cm^)) 
= (1.9 X 10-1,3.7 X 1028), (7.5 x 10-2,4.1 x 10^'^), and 
(1.1 X 10-1,8.0 X 102''), respectively. Here, dS is defined 
at the z = 0 and is parallel to the z-axis. Thus, the 
mass-to-flux ratio of the region normalized by its critical 
value (M/$)crit = (0.53/(37r))(5/G)i/2 is p = 21, 74, 
and 56, respectively. These values are much larger than 
the initial mass-to-flux ratio p = A. 


In figure m we show the evolution of the Ohmic, Hall, 
and ambipolar diffusion coefficients, po, Vh,, and pA, 
at the center of model Para as a function of the central 
density. The evolution of model Ortho and NoHall were 
almost the same. In pc < IQ-i'i g cm- 8 , pf^ is larger than 
PA and po and the gas rotation is significantly affected 
by the Hall current term in this region. The value of pu 
is much larger than the “critical va lue” for disk forma¬ 
tion ( thin black line) suggested by iKrasnopolskv et ai\ 
in this region. Although the pn decreases in 
Pc ^ 10-^8 g cm- 8 , the Ohmic and ambipolar diffusion 
alternatively play a role in pc > g cm- 8 , and the 

rotation is maintained in the high density region without 
magnetic braking. Note that IKrasnopolskv et a,l\ (120111 1 
only considered the Hall term and neglected other non¬ 
ideal effect. We remark that pA does not strongly depend 
on |B| around p g cm -8 although pA oc |B|2 in 

low density regions. 

Because of the conservation of angular momentum, 
the spin-up due to the Hall term at the center causes 
spin-down of the outer region, eventually, causing anti¬ 
rotation against the disk. In figure [SJ we show the cross 
section of the rotation velocity distribution of Model 
Para in the x — z plane at the same epoch of figure [TJ 
This figure clearly shows an anti-rotating envelope sur¬ 
rounding the forward rotating inner region. Since the 
anti-rotation of envelope is driven by torsional Alfven 
waves, the anti-rotating region expands with time and 
will propagate to the outside of the parental core. Thus, 
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Fig. 2.— Solid lines show the ratio of the sum of the centrifu¬ 
gal force and the pressure gradient force to the radial gravitational 

force, gi = | ^ I- Here, v^, p and $ are the rotation ve¬ 

locity, gas pressure and the gravitational potential, respectively. 
The dashed lines show the ratio of the centrifugal force to the ra¬ 
dial gravitational force, q 2 = \ \ ■ The dashed-dotted lines show 

g = 0.5. In the regions where the dashed lines are larger than the 
dashed-dotted lines, the gas is mainly supported by the centrifugal 
force. The top and bottom panels show the results of Models Para 
and Ortho,respectively. The epochs of each model are the same as 
those in figure [T] 

the angular momentum of the direction opposite to the 
disk would eventually be cast away to the interstellar 
medium. 


4. CONCLUSIONS AND DISCUSSION 

In this study, we investigated the effect of the Hall 
current term on the formation of circumstellar disk, by 
considering all non-ideal effects as well as the radiative 
transfer are considered. To our knowledge, this is the 
first study that simultaneously includes these physical 
processes in a three-dimensional simulation. 

We found that the disk evolution can be categorized 
into two cases depending whether the magnetic field and 
the rotation vector are parallel or anti-parallel. In the 
anti-parallel case, a relatively large (r > 10 AU) and 
massive disk forms simultaneously with the protostar for¬ 
mation; however, a disk with r > I AU does not form in 
parallel case. Thus, the parity of the magnetic field sig¬ 
nificantly changes the disk formation process which has 
not been paid much attention to so far. Since the parallel 



Fig. 3.— The time evolution of the mean specific angular mo¬ 
mentum of the inner region with p > 10“^^ g cm“^ as a function 
of the central density. The solid, dashed, and dotted lines show 
the results of Model Para, Ortho, and NoHall, respectively. 



Fig. 4.— Magnetic diffusion coefficients , 770, PH’, and pA at 
the center as a function of the central density in Model Para. The 
red line shows tjq , The green line shows rjn where the dashed line 
shows the region of rj^ < 0 and the solid line shows the region 
of rjH > 0, and the blue dashed-dotted line shows tja- The black 
thin line sho ws the “critical valu e ” of r in, 3 x 10^*^ Bc(cm^s“^) 
suggested by IKrasnopolskv et a g ii2oTiri above which the disk is 
formed in their simulations. Here, Be is the central magnetic field. 
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Fig. 5.— The cross section of in the x-z plane in Model Para. 
The epoch of the snapshot is the same as that in Fig. ^ 
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and anti-parallel properties are inherent and do not read¬ 
ily change in time, we expect that the spin-up (the mag¬ 
netic braking weakening) and spin-down (the magnetic 
braking strengthening) effects due to the Hall current 
term are also important in the subsequent disk evolution 
and the difference between the two cases is maintained 
or enhanced. Therefore, we suggest that the disk evolu¬ 
tion can be categorized into two cases. We may call the 
resultant parallel and anti-parallel systems as Ortho-disk 
and Para-disk, respectively. 

Our results predict that the bimodality in the disk size 
distribution spontaneously arises due to the Hall current 
term in typically magnetized molecular clouds. We tend 
to think that the disk size has an unimodal distribution 
according to the strength of the rotation and the mag¬ 
netic field of the cloud cores. However, as we have shown 
above, the Hall term changes the disk size according to 
the parallel or anti-parallel property of the cores. It is 
expected that about half of the molecular cloud cores 
have the parallel configuration and the others have the 
anti-parallel configuration because the Hall ter m would 
not p lay a role during the cloud core formation (|Wardlel 
12004) and there is no physical mechanism which distin¬ 
guishes the parallel and anti-parallel configurations. On 
the other hand, during the gravitational collapse of cloud 
core, the Hall term becomes effective and strengthens or 
weakens the magnetic braking. Therefore, the bimodal¬ 
ity of the disk size distribution spontaneously arises from 
the unimodal distributions of the rotation and the mag¬ 
netic field strength of the cores. 

An observational signature predicted from our results 
is the anti-rotating envelope in the Class 0 YSOs. At the 
end of the simulations, the anti-rotating envelope had a 
size of r > 100 AU and the rotation velocity of ~ 1 
km. Because the anti-rotation of the envelope stems from 
torsional Alfven waves, the anti-rotating region expands 
with time. We predict future observations will find an 
anti-rotating envelope against disk rotation in the Class 
0 YSOs. These observations will provide clear evidence 
that the Hall current term plays an important role for 
the evolution of circumstellar disks. 

In this paper, several simplifications were adopted and 
their influences should be investigated in future works. 
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strength of the Hall term, linearly depends on rjH- Thus, 
the simulations with the different models are necessary 
to confirm our results. The misalignment between the 
magnetic field and the rotation vector is another impor¬ 
tant issue. In our simulations, the initial rotation vec¬ 
tor and the magnetic field are in perfectly parallel or 
anti-parallel configurations. However, it is expected that 
they are mutually misaligned in the realistic c loud cores 
(|Hennebelle and Ciardil[200^ iJoos et alll2012() . Its effect 
on the disk formation with non-ideal effects should also 
be investigated. We used a rigidly rotating shell as the 
outer boundary condition. Because of the angular mo¬ 
mentum out-flux at the boundary, the total angular mo¬ 
mentum is non-conserved quantity in our simulations. At 
the end of the simulations, the total angular momentum 
of Models Para, Ortho, and NoHall within the boundary 
shell was 97, 94.3, and 95.5 % of the initial angular mo¬ 
mentum, respectively. Note that the difference of the an¬ 
gular momentum between models is expected because the 
Hall current term changes the angular momentum trans¬ 
fer rate near the boundary. Similar phenomena are also 
observed in the previous works. Previous simulations 
starting from the non-rotating core with the Hall term 
and outgoing boundary have a fini te angular momen¬ 
tum at the end of t he simulations (iKrasnopolskv et al\ 
[20nl: lLr^rZI[20Tl . We expect the treatment of the 
outer boundary condition would not change our results 
significantly because the crossing time of Alfven wave 
Across = R/va is larger than the free-fall time ts (in 
our initial condition, tcross/iff = 2.5) and the boundary 
mainly influences the relatively outer region within our 
simulation time tsim l.ltff. However, more sophisti¬ 
cated boundary conditions are desired. 

ACKNOWLEDGMENTS 

We thank K. Tomida and Y. Hori for providing EOS 
table. The computations were performed on the XC30 
system at CfCA of NAOJ. 


M. N. Machida, S.-i. Inutsuka, and T. Matsumoto, 

MNRAS 438, 2278 (2014), arXiv: 1307.1747 [astro-ph.SR] 

K. Tomida, S. Okuzumi, and M. N. Machida, 

ApJ 801, 117 (2015), arXiv:1501.04102 [astro-ph.SR] 

Y. Tsukamoto, K. Iwasaki, S. Okuzumi, M. N. Machida, and S.-i. 
Inutsuka, ArXiv e-prints (2015a), 

arXiv:1503.04901 [astro-ph.SR] 

K. Tomida, K. Tomisaka, T. Matsumoto, Y. Hori, S. Okuzumi, 

M. N. Machida, and K. Saigo, ApJ 763, 6 (2013), 
arXiv:1206.3567 [astro-ph.SR] 

M. Wardie and O. Ng, MNRAS 303, 239 (1999), 
astro-ph/9810468 

R. Krasnopolsky, Z.-Y. Li, and H. Shang, ApJ 733, 54 (2011), 
arXiv:1101.3018 [astro-ph.SR] 

Z. -Y. Li, R. Krasnopolsky, and H. Shang, ApJ 738, 180 (2011), 
arXiv:1106.2620 [astro-ph.GA] 

C. R. Braiding and M. Wardie, MNRAS 422, 261 (2012), 
arXiv:1109.1370 [astro-ph.SR] 

M. N. Machida and T. Matsumoto, MNRAS 413, 2767 (2011), 
arXiv:1008.0920 [astro-ph.SR] 




















6 


W. B. Dapp, S. Basu, and M. W. Kunz, A&A 541, A35 (2012), 
arXiv:1112.3801 [astro-ph.SR] 

K. Iwasaki and S. Inutsuka, MNRAS 418, 1668 (2011), 
arXiv:1106.3389 [astro-ph.GAJ 

K. Iwasaki and 8.-1. Inutsuka, in Numerical Modeling of Space 
Plasma Flows (ASTRONUM2012), Astronomical Society of 
the Pacific Conference Series, Vol. 474, edited by N. V. 
Pogorelov, E. Audit, and G. P. Zank (2013) p. 239. 

S. C. Whitehouse and M. R. Bate, MNRAS 353, 1078 (2004), 
arXiv:astro-ph/0406392 

S. G. Whitehouse, M. R. Bate, and J. J. Monaghan, 

MNRAS 364, 1367 (2005), arXiv:astro-ph/0509837 

Y. Tsukamoto, K. Iwasaki, and S.-i. Inutsuka, 

MNRAS 434, 2593 (2013a), arXiv:1305.4436 [astro-ph.SR] 

J. Wurster, U. Price, and B. Aylifie, MNRAS 444, 1104 (2014), 
arXiv: 1408.1807 [astro-ph.SR] 

V. Aiexiades, G. Amiez, and P.-A. Gremaud, Com. Num. Meth. 
Eng 12, 12 (1996). 

M. N. Machida, S.-I. Inutsuka, and T. Matsumoto, 

PASJ 63, 555 (2011a), arXiv: 1009.2140 [astro-ph.SR], 

T. Sano and J. M. Stone, ApJ 570, 314 (2002), astro-ph/0201179 
D. Semenov, T. Henning, G. Helling, M. llgner, and E. Sedlmayr, 

A&A 410, 611 (2003), arXiv:astro-ph/0308344 
J. W. Ferguson, U. R. Alexander, F. Allard, T. Barman, J. G. 
Bodnarik, P. H. Hauschildt, A. Heffner-Wong, and 
A. Tamanai, ApJ 623, 585 (2005), arXiv:astro-ph/0502045 

S. Okuzumi, ApJ 698, 1122 (2009), 
arXiv:0901.2886 [astro-ph.FPJ 

T. G. Mouschovias and L. Spitzer, Jr., ApJ 210, 326 (1976) 


T. Matsumoto and K. Tomisaka, ApJ 616, 266 (2004), 
arXiv:astro-ph/0408086 

M. N. Machida, S. Inutsuka, and T. Matsumoto, 

ApJ 670, 1198 (2007), arXiv:astro-ph/0702183 
C. Marois, B. Zuckerman, Q. M. Konopacky, B. Macintosh, and 
T. Barman, Nature 468, 1080 (2010), 
arXiv: 1011.4918 |astro-ph.FP| 

A. P. Boss, Science 276, 1836 (1997) 

S. Inutsuka, M. N. Machida, and T. Matsumoto, 

ApJ 718, L58 (2010), arXiv:0912.5439 [astro-ph.EP] 

M. N. Machida, S. Inutsuka, and T. Matsumoto, 

ApJ 729, 42 (2011b), arXiv:1101.1997 [astro-ph.SR] 

Y. Tsukamoto and M. N. Machida, MNRAS 416, 591 (2011) 

Y. Tsukamoto and M. N. Machida, MNRAS 428, 1321 (2013), 
arXiv:1210.0526 [astro-ph.SR] 

Y. Tsukamoto, M. N. Machida, and S. Inutsuka, 

MNRAS 436, 1667 (2013b) 

Y. Tsukamoto, S. Z. Takahashi, M. N. Machida, and S. Inutsuka, 
MNRAS 446, 1175 (2015b), arXiv: 1404.7271 [astro-ph.SR] 

M. Wardie, Ap&SS 292, 317 (2004), astro-ph/0307086 
M. Padovani, U. Galli, P. Hennebelle, B. Commergon, and 
M. Joos, A&A 571, A33 (2014), arXiv:1408.5901 
P. Hennebelle and A. Giardi, A&A 506, L29 (2009), 
arXiv:0909.3190 [astro-ph.GAJ 

M. Joos, P. Hennebelle, and A. Ciardi, ASzA 543, A128 (2012), 
arXiv:1203.1193 [astro-ph.SR] 


